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I. INTRODUCTION 

It is strongly believed nowadays that the universe is 
experiencing an accelerating expansion, and this is sup- 
ported by many cosmological observations, such as SNe 
la Q, WMAP ;2], SDSS [| and X-ray Q. These observa- 
tions suggest that the universe is dominated by dark en- 
ergy, which provides the dynamical mechanism for the ac- 
celerating expansion of the universe. Furthermore, they 
suggest that the dark energy equation-of-state parameter 
Wde might have crossed the phantom divide — 1 ^5|] from 
above in the near past. 

In order to provide theoretical paradigms for the de- 
scription of dark energy, one can either consider theories 
of modified gravity 0, or field models of dark energy. 
The field models that have been discussed widely in the 
literature consider a cosmological constant 0, a canon- 
ical scalar field (quintessence) [1], a phantom field, that 
is a scalar field with a negative sign of the kinetic term 
P, 10], or the combination of quintessence and phantom 
in a unified model named quintom (Tl] |. The advantage 
of quintom models is that they are capable of describing 
the crossing of the phantom divide, since in quintessence 
and phantom models wde can indeed vary, but it al- 
ways remains on the same side of the phantom bound 
(woE > — 1 in quintessence and woe < — 1 in phantom 
scenarios). 

On the other hand, dark energy models with nearly 
flat potentials have been shown to present interesting 
cosmological behavior, especially in the case where wde 
is around —1. Although there is not a concrete proof, 
there are many arguments indicating that nearly flat po- 
tentials are a natural way of acquiring wde ~ —1, if one 
desires to avoid sophisticated, and difficult to be justified, 
constructions. In the case of quintessence and thawing 
quintessence the corresponding conclusion has been ac- 
quired in various ways [12], while in the case of phantom 
cosmological paradigm in 13]. Finally, note that flat 
potentials, which keep the variations of the scalar fields 



from their initial to their present values small, can also be 
efficient to avoid unknown quantum gravity effects [3] ■ 

In [l^ the authors showed that all quintessence models 
with nearly flat potentials converge to a single function 
for w{a), with a the scale factor, which can be approxi- 
mately given analytically. Similarly, in [T6| the authors 
result to such a limiting behavior for the phantom sce- 
nario with nearly flat potentials. Both works share the 
general feature of quintessence and phantom paradigms, 
that is they cannot describe the —1 crossing, remaining 
either above (quintessence) or below (phantom). 

In the present work we are interested in investigat- 
ing the behavior of the combined case, that is we study 
quintom models with nearly flat potentials. Note that we 
do not assume complete dark energy domination, since 
matter is always non-negligible. We provide analytically 
an approximated universal behavior for wde{z), which 
can naturally describe the crossing of the phantom di- 
vide from above in the near past. This feature could 
bring the model at hand closer to observations. 

In section|lT]we construct general quintom models with 
nearly flat potentials, and we extract the expression for 
wde{z)- In section Hill we examine its behavior and we 
discuss the cosmological implications. Finally, in section 
IIVI we summarize our results. 



II. QUINTOM MODELS WITH NEARLY FLAT 
POTENTIALS 

Throughout the work we consider a flat Robertson- 
Walker metric: 

ds^ ^dt^ ~a^{t)dx^, (1) 
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with a the scale factor. The action of a universe consti- 
tuted of a canonical (f> and a phantom a fields is [T]| : 



d x\J—g 



(2) 



where we have set SttG = 1. The term Cm accounts for 
the matter content of the universe, which for simplicity 
is considered us dust. The Friedmann equations and the 
evolution equation for the canonical and phantom fields 
are [ul: 



= \ ^PM + P<t> + Pa] . 



(3) 
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^ ^3 . 



^ + 2p0 + 2pa 



V^{cl>) + V,{a) , (4) 



da 



0, 



(5) 



(6) 



where H = a/a is the Hubble parameter. In these ex- 
pressions, p0 and p0 are respectively the pressure and 
density of the canonical field, while and pa are the 
corresponding quantities for the phantom field. Finally, 
Pm is the density of the matter content of the universe. 

The energy density and pressure of the canonical and 
the phantom fields, are given by: 



and 



P0 
Pa = 

Pa = 



\4'^ - VM 



-CT^ + V,{a) 



(7) 



K(<t). 



(8) 



As usual, the dark energy of the universe is attributed to 
the scalar fields and it reads: 



1 



DE 



(9) 



Finally, the equation of state for the quintom dark energy 
is [11: 



WDE 



+Pa 
P4, + Pa 



2[V^(C^)+V^{CJ)] 



2[V^{(p)^V„{a)] 



(10) 



Using the definitions for the energy densities and pres- 
sures ([7]), ([8]), the Friedmann equations (l3|),([l]) can be re- 
written as: 



H = -\\PM 



PM + ^^' + v4<j))-h' + vA<j) 



(11) 

(12) 



Lastly, the equations close by considering the evolution 
of the matter density: 



Pm = ~iHp 



M- 



(13) 



In order to provide an analytical expression for wde '^^ 
have to transform the dynamical system (fSl) . ([6|) . (fTTj) . (fT2|) 
into an autonomous form flT'l. This will be achieved by 
introducing the auxiliary variables: 



V6H 



VV4^)+Vaia) 



y^(0) + K(a) y^TT"^ 



5 = 



' /V^W j ^ dV„(a) 



V^i<l>) + V4a) 



(14) 



(15) 



(16) 



(17) 



together with M = Ina. Thus, it is easy to see that for 

r dF 

dM- 



every quantity F we acquire F — Hj^. 



Using these variables, we result to the following sys- 
tem: 



dx 
dM 



^x{l + x'~y')-3x+Jhy' (18) 



dy 
dM 



'-y{l + x^-y^) 



5 xy 



(19) 



d\ _ Axiv,x\,xS) dV^(<j>) Bx{v,x\,xS) d'^V^{<p) 

dM v^i4>)+v„{a) d4, y^(0)+v„(<T) d</.2 -1- 

Cx{y,x\,xS) dV„(a) Dx{y,x\,x5) d^V„(a) 

V^(<t>)+V„{a) da V4,(4>)+V„(a) da'^ 



dS _ As(y,x\,xS) dV^i^) ^ Bs{y,x\,x5) d^VV(0) 



(20) 



dM — V^(4,) + V^(a) 



V4,(4>) + V^{a) ~~dp ^ 



C/i{y.x\,xS) dV^{a) . Ds{y,x\,xS) d Vcyja) /rj-i \ 

V^(4>)+V„(a) da + V^{4>) + V„{a) da^ ' ^'^^1 
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The complicated functions A\, B\,C\, D\,As, Bs,Cs, Ds 
can be straightforwardly calculated through differentia- 
tion of (fT6|) and (fT7|) . and elimination of 4> and & in terms 
of X, y, A, 5. However, their exact form is not needed for 
the purpose of this work. 

In terms of these auxiliary variables Qde from relation 
(O can be written as: 



DE 



y 



(22) 



Similarly, using (jlOp we can find the corresponding rela- 
tion for wde- However, it is more convenient to define a 
new variable Q = 1 + wde, which simply reads: 



C = 1 + Wde 
Thus, we can easily see that: 

2 C^DE 



2x^ 



(23) 



The first approximation is that \(\ ^ 1, i.e wde is close 
to —1. This assumption is justified by observations both 
at present time and in the recent cosmological past. The 
second assumption is that the potentials of the model 
are nearly flat, which is the case of interest of the present 
work. In particular we assume: 



< 1 

< 1 

< 1 

< 1. 
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dV^ict)) 




+ V„{a) 


del) 
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dV^ia) 




+ V,{a) 


da 




1 


d^V^{cp) 




+ V„{a) 


d(j? 




1 






+ V^{a) 


da^ 



(29) 



Therefore, the variables A and S can be assumed to be 
constant, and equal to their initial values: 



y 



DE 



1 



(24) 



Therefore, by differentiating these expressions with re- 
spect to M = Ina and using P^ . p^ . we acquire the 
autonomous equations for Qde and 



dflr 



dM 



ZVIde{1 - 0(1 - ^de) + VIdeVK^l 



2 



dC_ 

dM 



-3C(2 - C) + A(2 - O^K^DE - 



where 



= \-6. 



Finally, we obtain: 



_dQ_ 



dM 



-3C(2-C) + A(2-C)N/3CnaE-V3CnaE(l-|)C^ 

3J1l3i3(l-C)(l-fiDE)+^2£>E\/3a2DE(l-§)« 



(25) 



(26) 



(27) 



(28) 



A = Ao 



1/^(0) + K (a) 



Po- 



dV^{<t>) 
d4, 



(30) 



V^{<^) + V.{a) 



(31) 

with </)o, ctq the initial values and 0o, a^o the initial deriva- 
tives of the fields just before they begin the roll-down to 
the potential (equivalently the corresponding values in 
the near cosmological past). Thus, |Ao|,|<5o| <C 1, too. 
Note that at first sight, one could think that this ap- 
proximation would brake down if (jjQ and ao are equal, 
or more generally if (j) and a become equal at any time. 
However, this is not happening since in this case x in (jl4p 
and C in ([23]) also go to zero, with the hmits AoVC, (^o\/C 
not only regular but proportional to the first potential 
derivatives. That is, equation (j28p becomes even sim- 
pler. Finally, we mention that in the case where 4> and a 
become equal at some time, the derivatives -j^ and^^ 
remain small since they also depend on the combinations 
xX and xS. 

Keeping terms up to lowest order in ^, Aq, Sq, 
yields: 



Equation (|28|) . which gives the autonomous evolution 
for C (i.e for the equation-of-state parameter) is exact 
and takes into account the full dynamics of the system 
evolution, which is also determined exactly by equations 
p8| - (|2ip . In order to proceed to the extraction of ana- 
lytical solutions we have to make two assumptions. This 
will allow us to bypass the details of specific models and 
describe the general behavior of quintom scenarios with 
nearly flat potentials. Fortunately, as we are going to see 
in the next section, the error of our approximated solu- 
tion comparing to the numerical elaboration of the exact 
system is small. 



dC 



dn 



DE 



-6C + 2\o^3CnDE - ^3C^deC0o 
3f^D£;(l — ^de) ' 



(32) 



where ^0 = Xq — So- Equation (|32p can be transformed 
into a linear differential equation under the transforma- 
tion s = and can be solved exactly. The solution for 
C = l-f Wde is: 



1 -I- wde{^de) 



^jmo^ tanh 



2V3 



3+29oAo 



In 



(33) 
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Equation ([33)) . along with the corresponding result for 
w{z) derived below, is our main result. It shows that 
for sufhciently flat potentials, all quintom models with 
wde ~ — 1 at present time, approach a single generic be- 
havior. The specific role of the potentials is to determine 
the small constants Ao and Oq. Finally, we mention that 
(I33p gives always real values for wde as expected, due 
to the specific combination of the two terms that may be 
imaginary. 

In order to express our result in a form more suitable 
for comparison with observations we will use (|25p up to 
lowest order in ( and Oq, in order to acquire flnEiM), 



i.e rioEiz), since e 



M 



{1 + z) ^ (we set the present 



value of the scale factor ao = 1). In this approximation 
the solution is: 



(34) 



where floEo is the present-day value of Qde- Note that 
equation (j34p coincides with the expression for Qde in 
[18j . where the authors obtain in a different framework 
but under the assumption wde ~ ^ 1- 
Substituting ((34l) into ((33)) we obtain: 



1 + wde{z) = 



1 + (1 + zf 



1 



n 



- 1 



DEo 



-,1/2 



/ 3 + 29qXq 



/ 00 / 3 + 2goAo 



In 



1/2 



1/2 



1 



(35) 



We can expand ([55)1 with a — {1 + z)^-^ around a ~ 1, 
and acquire a linear relation of the form woEio) = 
WDEo + WDEai^ ~ '^); ii^ accordance with Chevallier- 
Polarski-Linder parametrization fl^. In this case the 
parameters wdeq ^ind woEa sre independent and are 
given as functions of Aq, Oq and CIdeq- This is an ad- 
vantage comparing to the simple quintessence [Tsj and 
simple phantom [16(| cases where such a procedure leads 
to a dependence of wde^ on wdEoj in contrast with the 
imposed parametrization. The inclusion of two fields, 
and thus of additional degrees of freedom, in the model 
at hand, restores the linear parametrization in its cor- 
rect form, that is with two independent parameters. Per- 
forming the aforementioned expansion imposing riu Eq = 
0.73 ± 0.03, with |Ao| < 0.5 and l^ol < 0.5, we obtain 
-1.08 < WDEo < -0-92 and -0.14 < wde^ < 0.14. 
These limits are narrower than those arising from ob- 
servations P, [2j H, 0| , which was expected since in our 
analysis we have assumed that WDE{a) is close to —1 

(ICI«i). 

Before closing this section let us make some comments 
on the expressions ([55)1 and ([55)1 . If we desire to obtain 
the simple canonical field, that is the case of quintessence 
models, we have to set the quantities that are relevant to 



the cr-field to zero. Thus, ((30)) and ((3T)) imply Ao = (5o G 
M and therefore 6*0 = 0. In this case p3p gives: 



\2 

l+WDE = 



1 



DE 



1 



DE 



1 In 



1 + ^/^DE 



which is just the result obtained in [If 



1 — \/^DE 

In addition, 



in this case by inserting the present value wueq one 
can solve for Ao, and substituting in ()35p he can obtain 
wde{z) with only wdEq s-nd ^dEq as parameters. Note 
that for the quintessence scenario, the variable x defined 
in (fT4)) is always real and thus ((23)) implies that is 



always larger than -1. This is just what is expected for a 
quintessence model. 

On the other hand, if we desire to obtain the simple 
phantom model, then we have to set the quantities rele- 
vant to the 0-field to zero. Therefore, i5o = — Aq G I and 
thus ^0 = 2A € I. In this case one finds exactly the same 
expression ((55)) . with Ag being negative, which is just the 
result obtained in [l6j . Finally, in this simple phantom 
case, one can also eliminate Aq and acquire wde{z) in 
terms of wdEq and VLueo- We mention that now x is 
purely imaginary (see relation ((H)) ) and thus ((^ leads 
to Wde always smaller than —1. Again, this is what is 
expected for a phantom model. 

In the quintom model at hand, that is in the case where 
both the canonical and the phantom fields are present, 
Ao, <5o and 6q are purely real or purely imaginary, depend- 
ing of which field is dominant, and the three variables be- 
long to the same set each time. We stress that A, 5 and 
6*, as well just suitable variables which allows us 

to transform the system to its autonomous form, and are 
not related to any observables. Therefore, their purely 
imaginary character is just a statement of the dominance 
of the phantom field, that is it acquires a robust and phys- 
ical content. This becomes obvious by the fact that all 
observables are real. Indeed, we can easily see that Q,de 
and in ((^ . ((^ . ((55)) and ((55)) are always real in 
the case of purely imaginary 5, 9 and x. The only ef- 
fect is that Wde is below the phantom divide. But such 
a transition is exactly the motive of the present work. 

Finally, we mention that when a crossing of the phan- 
tom divide takes place, that is when a; crosses zero, our 
system remains regular, and this was expected since even 
the naively acquired singular behavior is not related with 
the initial cosmological equations but only with the trans- 
formation we use to solve them analytically (one could 
resemble the case at hand with the distinction between 
true singularities and coordinate ones in general relativ- 
ity). In particular, in this case all the limits not only 
do exist but are much smaller than one, and the au- 
tonomous equations become even simpler. This can be 
also confirmed by the observation that when « cr the 
cosmological equations ([5])-((n]) become significantly sim- 
pler, and thus the corresponding autonomous system is 
simpler, too. Lastly, the regular behavior of our solution 
procedure is also confirmed by the numerical elaboration 
of the next section. 
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III. COSMOLOGICAL IMPLICATIONS 

Let us now investigate the cosmological implications of 
the acquired results. First of all, we desire to check the 
accuracy of our approximated analytical solution ([55)1 . 
In fig. [T] we perform such a comparison. We have se- 
lected two potential choices, namely V^{(j)) — cj)^, K-(o') — 
cr^ (dotted curves) and V^{(1>) = Kt(o') = cr~^ 

(dashed curves) and we have numerically found the exact 
wde{^de) behavior of the cosmological system, for three 
combinations of the parameters Ao and 9o (fixing suitably 
the values of 0o a-nd cto). From top to bottom the group 




0.4 0.5 0.6 0.7 0.8 0.9 

DE 

FIG. 1: (Color online) A comparison between the exact result 



for wde{^de) for V^{(j>) = (j) , 14 (ct) 



(dotted curves), 



the exact result for wde{^de) for V^{cj)) = 4>~ , Va-(cr) — 
(dashed curves), and our approximated analytical result 
for wde{0.de) provided by relation II33V (solid curves). The 
group of curves from top to bottom correspond to Aq — 0.7, 
00 = 0.7, to Ao = 0.4, 00 = 0.4 and to Ao = 0.1, Oo = 0.4. 
The error of our analytical result increases with the increase 
of \o and 6o, and with the increase of the distance of wde 
from —1. 



imation scheme brakes down. But these cosmological 
scenarios are beyond the purpose of the present work. 
Finally, note that independently of the subsequent cos- 
mological evolution, negative power-law potentials [2^ 
can fulfill the nearly-flat potential conditions ((29l) if (f>o is 
sufhciently large. This feature shows that any potential 
can give rise to the type of models discussed here, as long 
us they satisfy (HH), and this was also shown in [l^ for 
the case of thawing quintessence scenario. 

Having determined the applicability area of our ap- 
proximated scheme we proceed to specific cosmological 
scenarios. In fig. [5] we depict wde{^de), given by (|33p . 
for five different models. Firstly, with the solid and the 
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FIG. 2: (Color online) wue{S^de^ provided by relation S33\) 
for various cosmological paradigms. The solid and the dashed 
curves correspond to simple quintessence and simple phantom 
models, with Ao = 0.5 and Ao ~ —0.5 respectively. The dotted 
curve (quintomi) corresponds to Ao — 0.1, 6o — 0.4. The 
sort- dashed- dotted curve (quintom2) corresponds to Ac = 0.3, 
6o = 0.1. The dashed-dotted curve (quintomS) corresponds to 
Ao = 0.3i, Oq = O.li. 



of curves correspond to Aq = 0.7, Oq — 0.7, to Ao ~ 0.4, 
6*0 = 0.4 and to Aq = 0.1, = 0.4. In addition, for the 
same parameter values we have used relation p3p to ob- 
tain the approximated analytical behavior (solid curves). 
As can be seen, the deviation from the exact solution 
is larger for larger distance of wde from —1, which was 
expected since we have approximated |1 -I- wue\ to be 
much smaller than 1. Secondly, we observe that the error 
of our approximated result increases with the increase of 
the parameters Aq and i-e from bottom to top, which 
was also expected since we have assumed |Ao|, |6'o| <C 1. 
Thus, fig. [1] constraints the applicability range of our ap- 
proximated analytical solution to —1.1 < wde < —0.9 
and I Ao 1, 1^0 1 ^ 0.5. We mention that for higher param- 
eter values or for models with wde larger than 0, the 
deviation of our approximated solution from the exact 
cosmological evolution can be dramatic, and our approx- 



dashed curves we present the simple quintessence and 
the simple phantom scenarios, that is when ^o = and 
Aq = 0.5, Aq = —0.5 respectively. These results co- 
incide with those of [lB| and [l^, noting that in these 
works the authors draw wde{cl) instead of wde{z). As 
it was known, the simple quintessence and the simple 
phantom models cannot describe the transition through 
the phantom divide —1, and wde remains always on the 
same side of this bound {wde > —1 for quintessence 
and Woe < — 1 for the phantom scenario). However, the 
combined consideration of both models can indeed de- 
scribe the — 1-crossing. In fig. [5] we depict WDEi^DE) 
for three different quintom models. The dotted curve 
(quintomi) corresponds to Aq = 0.1, = 0.4, the sort- 
dashed-dotted curve (quintom2) corresponds to Aq = 0.3, 
00 = 0.1 and the dashed-dotted curve (quintom3) corre- 
sponds to Ao = 0.3i, 9o = O.li. As we observe, all three 
models present the phantom-divide crossing. 
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In order to acquire a more transparent picture, in 
fig. [3] we depict W£,e{z) (given by relation ([55]) with 
^DEo ~ 0.73) for the same cosmological models of fig. [2] 
As we observe, although all three examined quintom 
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dominance of one of the fields is sufficiently strong, then 
wde{z), although moving towards —1, it will never suc- 
ceed to cross it. 

Let us close this section with a quantitative discussion. 
Current observations of wde{z) [i|, 0, H, 0] can accept 
the description of quintom models with nearly flat po- 
tentials presented here. The most significant advantage 
of this scenario is the capability of describing the cross- 
ing through the phantom divide from above to below, 
in the near cosmological past, if such a crossing will be 
retained by future and more exact observations. In ad- 
dition, even without a crossing, the model at hand can 
describe the t(;D£;(z)-evolution in close agreement with 
current observational limits. However, we mention that 
such a description is quantitatively trustworthy if the cur- 
rent value wdEo is larger than —1.1 and smaller than 
—0.9. On the contrary, Hwdeo ^5 —1-1 or wdEq ^ —0.9, 
then the errors of our approximated analytical solution 
(|55)l become relatively large and the results have to be 
considered only qualitatively. 



FIG. 3: (Color online) woEiz) provided by relation \35\l with 
^deq ~ 0.73, for various cosmological paradigms. The solid 
and the dashed curves correspond to simple quintessence and 
simple phantom models, with Aq ~ 0.5 and Aq ~ —0.5 respec- 
tively. The dotted curve (qmntoml) corresponds to Ao = 0.1, 
6o = 0.4. The sort-dashed-dotted curve (quintom2) corre- 
sponds to Ao = 0.3, Oq = 0.1. The dashed-dotted curve (quin- 
tomS) corresponds to Ao — 0.3i, 9o = O.li. 



models can describe the phantom-divide crossing, only 
quirnoml and quintom2 present the desired behavior, 
that is wde{z) crossing —1 from above, resulting to 
wdEo < — 1 at present, as it might be the case accord- 
ing to observations. We remind that according to pO|) 
and (|3ip . imaginary values correspond to initial phan- 
tom dominance, while real parameter values correspond 
to initial canonical- field dominance. Thus our framework 
reveals a way of determining the specific features of a 
quintom model. In particular, the crossing through — 1 
requires a sign change of x'^ (see relation ([12])), and thus 
a change in which field's kinetic term is larger (see defini- 
tion (HH)). If initially (with "initial" having the meaning 
of the beginning of slow-roll as considered in [l^, [l^) is 
the canonical field that is dominant and at some point the 
phantom field dominates the evolution, then the crossing 
of —1 is from above to below. This case corresponds to 
quintomi and quintom2 paradigms of fig.O On the other 
hand, if initially the evolution is dominated by the phan- 
tom field, and progressively by the canonical field, then 
the crossing through —1 takes place from below to above. 
This is the case of quintomS scenario of fig. [3] There- 
fore, in order to describe the possible crossing of Wde{z) 
one has to consider an initial canonical field dominance, 
with a subsequent increase and dominance of the phan- 
tom field. These conditions are necessary for a certain 
— 1-crossing, bur are not efficient. That is, if the initial 



IV. CONCLUSIONS 

In this work we examine quintom models with nearly 
flat potentials, without the assumption of complete dark 
energy domination. In the case where the dark energy 
equation-of-state parameter wde close to —1, we pro- 
vide analytically an approximated universal expression 
for wde{z) for all such models. This expression depends 
only on the initial conditions (beginning of slow roll), i.e 
on the values of the potentials and their derivatives at a 
specific point, and on the values of the field kinetic terms 
at the same point. This feature arises because, due to 
the potential flatness, the fields never roll very far along 
the potentials in order to "feel" the rest of their shape. 

Contrary to the case of simple quintessence or simple 
phantom models, where wuEiz) is always on the same 
side of the phantom divide (wde > — 1 for quintessence 
and Wde < — 1 for the phantom scenario) the quintom 
paradigm allows for a description of the transition 
through —1. In addition, we provide the necessary 
conditions for a specific such crossing. In particular, if 
initially the universe is dominated by the canonical field 
then the subsequent evolution can bring the phantom 
field domination and thus the —1-crossing from above 
to below. On the other hand, if initially is the phantom 
field that dominates then the evolution can lead to 
— 1-crossing from below to above. Thus, a not- very- 
strong initial dominance of the canonical field can lead 
to a cosmological evolution where wde{z) crosses — 1 
from above to below in the near past. In conclusion, 
the determination of the impact of the two fields, can 
provide a WDE{z)-eYo\ntion in agreement with current 
observational limits. 
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